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A CLASS OF PERIODIC ORBITS OF SUPERIOR PLANETS* 

BY 

F. R. MOULTON 
1. Irvtroduelion. 

The periodic solutions of the problem of three bodies fall into many classes 
requiring somewhat different modes of treatment. One of these classes f is 
that in which a small, or infinitesimal, body moves near another of finite mass 
and is subject to the disturbing influence of a distant large mass which moves in 
a circle about the finite mass except in so far as it is disturbed by the small body. 
Moreover, the orbit of the small body is circular except as it is disturbed by 
the distant large mass. 

In this paper the problem of a distant infinitesimal body moving subject to 
the attraction of two finite masses relatively near each other and revolving about 
their common center of gravity in circles is considered, and the existence of 
certain periodic solutions is established. As the finite bodies approach coinci- 
dence the orbit of the distant infinitesimal body approaches a circle. The re- 
striction that the distant mass shall be infinitesimal is no more necessary than 
it was in the earlier problem, but it is made for the sake of simplicity. 

The mode of treatment here is similar to that adopted in the earlier investi- 
gation. For the value zero of the parameter n the problem reduces to that of 
two bodies which admits a circular orbit as a periodic solution. The existence 
of the analytic continuation of this orbit as a periodic solution is proved, and 
direct methods of constructing the series are developed. It is shown how the 
integral can be used as a check on the computations, or to replace one of the 
differential equations in the construction of the solutions. 

The earlier paper had direct and important applications to the Lunar Theory 
and to the satellite theory in general. The results of this paper do not find 
applications in the solar system, except possibly in the case of the most remote 
dmall satellites of Jupiter and Saturn which are subject to the attractions of the 
central planets and the larger interior satellites, but they have some use in 
certain triple star systems. Their chief value is certainly that they cover a part 

* Read before the American Mathematical Society, Chicago, April 29, 1911. This paper 
was written while the author was a Research Associate of the Carnegie Institution of 
Washington. 

t Treated by the writer in tnese Transactions, vol. 7 (1906), pp. 537-577. 
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of the field of periodic orbits of the restricted problem of three bodies, and an 
examination of the whole field is necessary for the discussion of the evolution 
of these orbits with a continuous change of the parameters upon which they 
depend, for example the constant of the Jacobian integral which Dr. Hill and 
Sir George Darwin adopted in their important researches.* The results of 
such an analysis as is given here supplement the laborious numerical investiga- 
tions of Hill and Darwin. 

2. The differential eqvptions. 

Let mi and TO2 represent the masses of the finite bodies and F the gravitational 
constant. Suppose the infinitesimal body is projected in the plane of motion 
of the finite bodies; its orbit will be a plane curve. Let the origin of coordinates 
be the center of gravity of the system and the a^-plane the plane of motion. 
Let the coordinates of mi, 7re 2 , and the infinitesimal body be xi, y\; x it y 2 ; and 
x , y respectively. Then the differential equations of motions for the infinitesimal 
body are 



(1) 



<Px^_dU <Py dU 

dfi ~ dx ' dP ~ dy ' 

TT — ^ mi 1 & m i 



n= •(*-*!)»+ (y-yi)», r 2 =^(x-x t )*+(y-y 2 )\ 



Let 

(2) 



r=Vx*+tf, R=V / (x 1 -x 2 ) i +(y 1 -y 2 y, 



m ^ t> n ,/-*-, 1 m l 



Ri=Vx\ + y\ = -^—R, R 2 =Vxi±yl = --^-R. 
m\-\- m.2 mi+wi2 



Then equations (1) become in polar coordinates 

(dvV_ 
\dt) ~ 



cPr_ (dvV_dU 
d# r \dt) ~ dr 
(3) 

eP» dr<fo_ 1 dU 
r ja 1 * Jt ji ~ 



df^dt dt r dv ' 
We now develop the potential function U. From (1) and (2) we find 



* Hill, American Journal of Mathematics, vol. 1 (1878),.pp. 5-26, 129-147, 245- 
260; Darwin, Acta Mathematica, vol. 21 (1897), pp. 99-242; Mathematische 
Annalen, vol. 51 (1898-9), pp. 523-583. 
Trans. Am. Hath. Soc. 7 
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1-i 



U=—\l- — cos (t»-t»i) + {-) j 

+ — -[1 + — COS (»-Di) + ^y J J 
I) 

fc 2 (mi+m 2 ) . ^mima fl 2 /, 
= h -jU l + 3cos2(t>-«i)] 

p D 1 

H — ^o — -"[3 cos (t> — »i) + 5 cos 3(» — »i)]+ •• • j • 

Then equations (3) become 

/dv\ 2 . P(mi + m 2 ) ft 2 mim 2 ft 8 J ar , 

H — ^-5 [3 cos (» — »i) + 5 cos3(» — «i)]+ • • • | , 



ffr 

dt 2 



(5) 



d?v dr dv _ Wmim^ R 2 



|| sin 2(» — »0 



-\ ^ [sin(»-»i) + 5sin3(i)-»i)]H J. 

If the orbits of mi and m 3 are circles, as is assumed to be the case, equations 
(1) admit the Jacobian integral * 

(t)'+(*)W't-'J]- w - c ' 

where 

(6) iti = pi > = constant of integration. 

It follows that ni is the mean motion of the finite bodies and that »i = Wi ( t — < ) • 
The integral becomes in polar coordinates 

When the right members are neglected equations (5) have the particular 
solution 

/ON k l/m l + m 2 /■ , ^ / 4 4 \ 

(8) r=a, v= -, (< — to) = n( <— <<>)> 

where n is the angular velocity of the infinitesimal body in its orbit and U is 
an arbitrary constant. It will be supposed that n is given by the observations, 
* See Moulton's Introduction to Celestial Mechanics, p. 202. 
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or that its value is assumed, and that o is determined by the second equation 
of (8). The constant a has three determinations, only one of which is real. 
New variables p , 6 , and t and new constants p and M will be introduced by 
the equations 

r = a(l+p), v=n(t — to) + d, 

TV 

(9) — =n, (n — ni) (t - t ) = t , 

«i 

mi m 2 

= M. 



[m l + m 2 ] 2 
It follows from (6), (8), and (9) that 

(10) -=M 5 - 

a 

Then equations (5) become 

1p_ m , J m del* m 2 1 

dr 2 Ui_p; Ll-M d<J + (l-M) 2 d + p) 2 

= ( 1 _; ) y +p) ^ [ 1 + 3co S 2(r + g ] 

< u > + r 1 ;m 2 i 2(i^p) t 3cos(T+ ^ +5c ° s3(T+ ^ + -i 

m 2 — mi 3u' . . . . . „, v , i 

+ ^r^T sTiTl sin < T + e + 5 sin 3 r + 6 )]+ • • • I. 
mi + m 2 8(1 + p) J 

These equations are valid for the determination of the motion of the infinitesimal 
so long as | n | < 1 . The right member of the first one involves only cosines, 
and of the second one only sines, of integral multiples of ( t + 6 ) . The parts in 
the brackets proceed according to powers of /J , the coefficients of even powers 
of i^ in the first and second equations being cosines and sines respectively of 
even multiples of ( t -(- ) , and the coefficients of odd powers of ju 3 being cosines 
and sines of odd multiples of ( t + 6 ) . 

3. Existence of periodic solutions. 

Suppose p — /3, dpjdt =0,6=0, dd /dr = y at t = and let the solution 

of (11) be written in the form 

P=/(/3,7;r), 
(12) 

6 = <f>(0,y, T ). 

" = <PKP,y\r). 
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Now make the transformation 

(13) p= p', 0=-8', t=-t'. 

The resulting equations have precisely the form (11). Consequently their 
solutions with the initial conditions p' = 0, dp' Jdr = 0, 0' = 0, dd 1 jdr = y are 

p' = f(P,y, t')=/(0, 7 ;-t) = p , 

(14) 

e' = <f>(P, r, T') = <p(p,y,-T)= - e. 

Therefore, with the initial conditions adopted, p is an even function of t and 9 
is an odd function of r . The orbit is symmetrical both geometrically and in t 
with respect to an axis passing through the origin and rotating with any con- 
stant velocity. In particular, it is symmetrical with respect to the line joining 
the finite bodies. Such an orbit will be called symmetrical whether it is periodic 
or not. 

Now consider the conditions for a closed symmetrical orbit. Since the right 
members of (11) involve sines and cosines of integral multiples of r, sufficient 
conditions that p and in symmetrical orbits shall be periodic with the period 
2jic , j an integer, are 

d P !dr = f'(p, r ,JT) = 0, 
(15) 

6 = <p(P, y,j*) = 0. 

These conditions are necessary provided they are distinct. 

In order to examine the solutions of (15) it is convenient to use other para- 
meters than j3 and y . Suppose that at t = 

r=a(l+p) = a(l + «)(l-e), 

dr dp 
(16) dr =a -dl =0 ' *-°' 



dd p. VI + e 



1-M dt l_ M (l + a )l(l_«)l- 

It follows that a ( 1 + a ) and e are the major semi-axis and the eccentricity of 
the elliptic orbit which we should obtain if the right members of equations (11) 
were zero. Because of the well-known properties of the solutions in terms of 
these elements, the properties of the general solutions so far as they do not depend 
upon the right members of (11) are known. These properties will be important 
in solving the conditions for periodic solutions. 

Equations (11) are regular in the vicinity of m = 0, p=0, dp/dt = 0, 
= , dd jdr = for all t . It follows from Cauchy's general existence theorem 
on the solutions of differential equations that the solutions are analytic functions 
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of a , e , and /t 1 which are regular in the vicinity of a = e = /i* = 0; and it 
follows from PoincarS's extension* of Cauehy's results that the moduli of a, 
e , and /** can be taken so small that the solutions are regular while r runs through 
any finite preassigned range of values. We shall choose as the interval for r 
the range = r = 2jv and integrate (11) as power series ma, e, and /»* > van- 
ishing with a = e = /** = . That is, our results are the analytic continu- 
ation with respect to these parameters of the particular solution r = a, v= nt 
which we have when p = . The results may be written in the form 

P = Pi(«' e »/* i ;T), 

dp i \ \ 

-fc = Pj(«,c,/»*;t), 

(17) 

6= p 3 (ct,e,fi*;T), 

de ( i ^ 

^= p 4 (a,e,M s ;T). 

where pi , • ■ • , p< are power series in a , e , /** with t entering in the coefficients. 
The conditions (15) for a symmetric periodic solution become 

P2(a,e,n i ;jir) = 0, 
(18) 

pi(a,e,n*;jic) = 0. 

It will be shown that these equations can be solved for a and e as power series 
in /x* , vanishing with m* — , and converging for the modulus of m* sufficiently 
small. 

Since the right members of (11) carry a» v as a factor, the part of the solution 
depending on the right members will be divisible by m v • If the right members 
of (11) were zero and the solution were formed with the initial conditions. (16), 
the mean angular motion of the infinitesimal body in its orbit would be 

(19) 



(1-M.) (l+«) r 



Consequently it follows from the solution of the two-body problem that equations 
(18) have the form 

P2OV) = ev { sin vj* + e sin 2vjv + • • •} + ft*q i (a,e,n i ; jr) = 0, 
(20) r 1 

PsOV) = j YZT~J* ~ V 3 T ~ 2e sm V 3 T ~ '"] +M V ?s(a,e,M*;>) = 0, 

where the unwritten parts in the { } are sines of multiples of vjv and carry e* 
as a factor. 



♦Acta Mathematica, vol. 13 (1890), p. Id. 
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On referring to (19) it is seen that the first equation of (20) is divisible by n 2 and 
the second by n . After division by these factors the equations are still satisfied 
hya=e=n=0. Moreover the determinant of their linear terms in a and e 
is 

0, jr 



(21) A = 



! = - %?** 



§ jw , — 2;V j 



&V + 0. 



Therefore, besides the solution fi = , equations (20) have a unique solution 
for a and e as power series in ju* , vanishing with ^ = , which converge for the 
modulus of n* sufficiently small. These power series carry ^ as a factor and 
can be written 

(22) a = /i*iW), e = M *P 2 ( M *). 

On substituting these series in the right members of (17), which vanish with 
a = e = n* = , we have 

(23) P =M s <?i(M';r), d = ^QA^;r). 

The series Qi and Q 2 are periodic in t with the period 2jir because the conditions 
that the solutions shall have this period have been satisfied. Since (17) converge 
for all = t = 2jw if the moduli of a , e , and n h are sufficiently small, and since 
the expressions for a and e given in (22) vanish for n = 0, it follows that the 
modulus of n* can be taken so small that the series (23) converge for all t in the 
interval; and since they are periodic with the period 2jw, the convergence holds 
for all finite values of t . 

The integer j has so far been undetermined. When j is unity the periodic 
solutions exist uniquely and their period is 27r . When j is greater than unity 
the periodic solutions also exist uniquely. Since the periodic orbits for j greater 
than unity include those for j equal to unity, and since in both cases there is 
precisely one periodic orbit for a given value of fi* , it follows that all the sym- 
metrical periodic orbits of the class under consideration hare the period 2w . 

It follows from (0) and (9) that t + = v — n . Since in the periodic solution 
6 is periodic with the period 2x , the period of the solution is the synodic period 
of the three bodies. Hence, if we refer the motion of the infinitesimal body 
to a set of axes having their origin at the center of gravity of the system and 
rotating in the direction of motion of the finite bodies at the angular rate at 
which they move, and if the .r-axis passes through the finite bodies, the 
periodic orbit of the infinitesimal body, which has been proved to exist, will be 
symmetrical with respect to the .r-axis. Since by hypothesis a > R it follows 
from (6) and (S) that n x > n . Therefore even if the motion of the infinitesimal 
body is forward with respect to fixed axes it is retrograde with respect to the 
rotating axes. 
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It is supposed that the period of the finite bodies, and therefore n\ , is given 
in advance and remains fixed. The variation of the parameter ju* corresponds 
to a variation of the period of the infinitesimal body defined by n . If the motion 
with respect to fixed axes is forward n has the same sign as n\ and p i has three 
values, one of which is real and positive while the other two are complex. If 
the motion is retrograde /t J has three different values, one of which is real and 
negative while the other two are complex. Therefore for a given period there 
are six symmetrical orbits, three direct and three retrograde, and for small | p \ 
one direct orbit is real and one retrograde orbit is real, while in the others the 
coordinates are complex. This means, of course, that the corresponding solu- 
tions do not exist in the physical problem. The coordinates of the complex 
orbits are conjugate in pairs. For a certain value of \£ they may become equal, 
and therefore real, and for larger values of /** real and distinct, as Darwin found 
by numerical computations to be true for satellites in a special case in A c t a 
Mathematica, vol. 21 (1897), pp. 99-242. 

Transforming the integral (7) by (9), we find 

^'Y-un-u s J m dej (i+p) 2 f m del y? 

\dr) +Ui " p; ll-p dr] l-/x Ll-M drj-(l-„)»(l+p) 

2M a? f , 

( 24 ) +( i_ M) 2( 1+p) a (i[l + 3cos2(r + e)] 

+ nr«M M l J 3cos(T+g) + 5cos3(r+g)]+ • • • } - d, 
nil + 7W2 o (. 1 -+- p ) J 

where C = n\{\ — p) 2 a 2 C\. It follows from this equation and (23) that C\ 
can be expanded as a power series in ju* , vanishing with p} = . The term of 
lowest degree in p} after substituting (23) is pi . Therefore m* can be expanded 
as a power series in C\ . For C L = the three branches of the function are the 
same. Since a = R J p i the relation between C and C\ is 

(25) C = j = — p 1 1 + power series in p J . 

From this it follows that C = oo for m* = . Therefore the periodic orbits 
branch at C = oo , and there are two cycles of three each. 

4. Practical construction of the periodic solutions. 

It has been proved that the symmetrical periodic solutions under discussion 
are expressible in the form 

(26) P=llpi^, 0=f><M* 3 , 

i-4 <=4 
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where the p< and 0< are functions of t . Since these series are periodic and con- 
verge for all | jt** | sufficiently small, it follows that each p< and 0,- is separately 
periodic; that is, 

(27) Pi(t+20- P i(t), 0i(r+2ir)^e i (r). 

In every closed orbit there are points at which dp / dr = . We shall suppose 
that <o of (9) is so determined that this condition is satisfied at t = . It follows 
from this and the convergence of (9) for all | m* | sufficiently small that 

(28) ^ = = 4, •••,«,). 

In the symmetrical periodic orbits whose existence has been proved above we 
have also — at t = . But we shall not now impose this condition because 
the general periodic orbits include those which are symmetrical; and if in the 
construction it appears that the conditions for symmetry are a consequence of 
those for periodicity, we shall have proved that all the periodic orbits of the class 
under consideration are symmetrical. Such will be seen to be the case. 

Equations (26) are to be substituted in (11), arranged as power series in ju J , 
and the coefficients of the several powers of ^ set equal to zero. The coefficients 
of /«* set equal to zero give the equations 

(29) -j^-O, -^ = 0. 

The solutions of these equations satisfying (27) and (28) are 

(30) Pi = a*> i = b i , 

where o« and b t are constants which are so far undetermined. 

The coefficients of m s , • • • > M 1 are the same as (29) except for their sub- 
scripts, and their solutions satisfying (27) and (28) are similarly 

(31) Pi=a it 0, = bj (i-5, ••-,»), 

where the oj and 6,- are all undetermined constants. 
The coefficients of m ¥ give the equations 

^- =3« 4 -^[l + 3cos2r], 

(32) (Pfc, ZM . _ 
-ST- — J" ■»*• 

In order that the solution of the first of these equations shall be periodic we 

must impose the condition 

M 

(33) a t = -j , 
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which uniquely determines the constant o 4 . Then the solutions of (32) satis- 
fying (27) and (28) are 

Pio = aio + •& M cos 2t , 
(34) 

0io = feio + %M sin 2t , 

oio and feio remaining undetermined. 
The coefficients of /tt* are 

from which it follows that 

(35) Of, =. 0, pu = on, On = feu, 

where on and fen are as yet undetermined. 
The coefficients of it* are 

d?Pi2 , M(mi- mjt) fo 

-d? = 3fl<+ 2(mi+ms) [3cosr+5cos3r], 

<ffli2 3 Jkf (mi — rm) . . 

~d?- = S mi + m, [smr+5sin3r]. 

Imposing conditions (27) and integrating, we find 

O6=0, 

M(mi— rrii) . 

(36) m = aii "2(mi + m 2 ) 13 C ° S T + * ^ 3t] ' 

, ZM ( mi — mj ) , . , K . „ , 

«u = 612 fiT— , - x [sin t + $ sin 3r] , 

o v. mi -J- t»2 ; 

where a w and feu remain undetermined. 
In a similar way we find from the coefficients of j» v 

(37) 07 = , pu = ois + fikf cos 2t , 0i» = fei 3 + &M sin 2t , 

an and feu remaining undetermined. 

So far all the fey have remained arbitrary and it is necessary to carry the in- 
tegration one step further in order to see how they are determined. The coef- 
ficients of /x v are 

-jjr = 3a 8 - Za\ + $M [36, sin 2t + 2a, + 6a 4 cos 2t1 , 
(38) 

#014 #010 „.„, - o • O 1 

-j-j- = — a« -t] — 3Af [ 64 cos 2t — 2a« sm 2t] . 
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Substituting the values of a 4 and 0io from (33) and (34), imposing the con- 
ditions (27) and integrating, we find 



3M 2 
as - - 16 , 

(39) Pl4 = a u - |&4 M sin 2t - T 9 ff iP cos 2r , 
0u = b u + \b t M cos 2t - £|M 2 sin 2t . 

Imposing the condition (28) for j = 14 , we have 

64 = 0, 

(40) Pu = a u — ^M 2 cos 2t, 

e u = b u - §§M 2 sin 2t. 

From the coefficients of /» v it is found in a similar way that 
h= 0, 

Pl5 = Ol6 - f -M — , ~ I COS T + 2 B T COS 3t ] , 

On = b n + UI ^'7^ [sin r + 2 5 T sin 3r] , 



(41) 



au and 615 remaining arbitrary. 

It will be observed that so far as the computation has been carried the coef- 
ficients of the pj are cosines of integral multiples of r and the coefficients of the 
0,- , except for the undetermined additive constants, are sines of integral multiples 
of t. In the computation of p, the periodicity conditions have uniquely de- 
termined cij-t , and the condition dpj jdr = at t = has required that £>/_io = . 
It will now be shown that these properties are general. Suppose p 4 , • • • , p n , 
$4, • • ■ , 0» have been computed and that the coefficients are all known except 
a„_6 . • • • > a« which enter additively in p»_6 , • • • , p n respectively, and &*_9 , 
• • • , b n which enter additively in 0„_9 , • • • , 0» respectively. For the deter- 
mination of Pn+i and 0„ + i we have 

"^T = 3a^_ B + flf 6-. sin 2r + F n+1 ( r ) , 
(42) 

where F n+ i ( r ) and G n +i ( t ) are entirely known functions of r . It follows 
from the assumptions respecting p 4 , • • • , p„ , 4 , • • • , 0* and the properties of 
equation (11) that F„+i ( t ) is a sum of cosines of integral multiples of t and that 
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G, l+ i ( t ) is a sum of sines of integral multiples of t . Hence we may write 

F B+1 (r) = S4" +1) cosir, 
G B+1 (r) = 2fl? +1) sin jr. 

In order that the solution of the first equation of (42) shall be periodic we 
must impose the condition 

(43) 3a n _ 5 +^; +1) = 0, 

which uniquely determines a„_s. 
After equation (43) has been satisfied the solution of the first equation of 

(42) is 

Pn+i = a„+i - |lf6 n _ 9 sin 2r + 2a? +1) cos jr , 

(44) ! 

a J ~ p^j 

The condition dp / dr = at t = makes it necessary to impose the condition 

(45) 6,-9 = 0. 

Then p„+i is completely determined except for the additive constant a n +i, and 
is a sum of cosines of integral multiples of r . 
The solution of the second equation of (42) is 

0,+i = 6„+i 4- 2#" +1) sin jr, 

(46) ! 

fl(»+l> _ _ _ RO+U 
PJ j2 °i 

Hence ^ n+ i is a sum of sines of integral multiples of t, except for the unde- 
termined constant 6 n+ i which must be put equal to zero in order to satisfy the 
conditions on p n+n . These results lead, by induction, to the conclusion that 
the Pi and jf j=4, • • •, oo, are sums of cosines and sines respectively of integral 
multiples of r Whose coefficients are uniquely determined. 

It follows from the properties of the solutions which have just been established 
that not only is dp / dr — at r = but also 6(0) = . Therefore these 
periodic orbits are the symmetrial orbits whose existence was established in § 3. 
In the construction it was not assumed that the orbits were symmetrical, and 
since this property was a necessary consequence of the periodicity conditions 
it follows that all periodic solutions which are expansible as power series in /** 
are symmetrical. It is easily shown by direct consideration of the construction 
of periodic solutions that they cannot be expanded as power series in ju Itf except 
when j is a multiple of 3 , and that then they reduce to those found above. 
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5. Application of the integral. 

The differential equations admit the integral (24) which, for brevity, can be 
written in the form 



'("&•• S".*)-«- 



It follows from the form of (24) and the expansion (26) that the left member 
of this equation can be developed as a power series in yl , giving 

(47) F = F + F u S + F 2l j*~+ ■-■+F n -J+--- = 0. 

Since the p,- and 6j are sums of cosines and sines respectively of integral 
multiples of t , and since dp / dr enters in (24) only in the second degree and 9 
only in even degrees, it follows that the F,- are sums of cosines of integral multiples 
of t . Equation (47) is an identity in ^ , whence 

F„ = xOf cos jr = (n = 0, ...,oo). 

Since these equations hold for all values of t , we have 

(48) Of=0 (n,j = 0,...,oo). 

The Of are functions of the elf, • • • , elf and /3/ 0) , ••• ,0f. Hence equations 

(48) can be used as check formulas on the computation of the coefficients of the 
solutions. 

Equations (48) can be used in place of the second equation of (11) for the 
determination of the (iff, the coefficients of the trigonometric terms in the 
expression for 0« . Suppose pt , • • • , p»_i and 04 , • • - , 0»-i have been deter- 
mined except for additive constants in pn-t, • • • , p„_i. It follows from (24) 
that F n is 

(49) '-£+*■(•■£•*■£) <"* -»■ 

where P n is a polynomial in the arguments indicated. Consequently equations 
(48) are of the form 

(50) Of=(ff+Vf(<$\ff) = (,-4,. ..,„_!), 

where the Bf are known polynomials in the elf and the fff . Therefore equa- 
tions (50) uniquely determine the fff 

The University op Chicago, 
April 29, 1911. 



